1^2 Topological Insulators in Ultracold Atomic Gases 
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We describe how optical dressing can be used to generate bandstructures for ultracold atoms 
with non-trivial Z2 topological order. Time reversal symmetry is preserved by simple conditions on 
the optical fields. We first show how to construct optical lattices that give rise to Z2 topological 
insulators in two dimensions. We then describe a general method for the construction of three- 
dimensional Z2 topological insulators. A central feature of our approach is a new way to understand 
Z2 topological insulators starting from the nearly-free electron limit. 
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Recently it has been realized that the electronic struc- 
ture of materials is much richer than had previously been 
thought. The apparently simple band insulator admits 
an intricate topological classification, and can behave as 
a "topological insulator" which is gapped in the bulk but 
carries gapless, metallic, states on its surface [l|, The 
classic example is the integer quantum Hall effect (IQHE) 
of a two-dimensional (2D) electron system in a magnetic 
field. In this setting, of 2D with broken time-reversal 
symmetry (TRS), band insulators are characterized by 
the integer- valued Chern number [3], which determines 
the number of metallic surface states. Presently, interest 
is focused on settings that preserve TRS, since unlike the 
IQHE, topological insulators are possible in both 2D and 
3D. They are classified by a Z2 invariant : band in- 
sulators either are trivial, or are non-trivial and support 
a metallic surface state with interesting spin-structure. 

While the effects of interactions are rather well un- 
derstood in the IQHE case, our understanding of strong 
correlations in Z2 topological insulators, especially in 3D 
is yet in its infancy. Ultracold atomic gases might hold 
much promise in this direction: they allow the effects of 
interactions and other perturbations to be explored in a 
controlled way0, S]- To take advantage of this in the 
topological insulator context, it is a crucial prerequisite 
to find a way to simulate bandstructures with nontriv- 
ial Z2 invariant. The purpose of this Letter is to solve 
this problem in both 2D and 3D in a simple, generic and 
experimentally feasible fashion. 

The standard technique for imposing bandstructures 
on ultracold atom gases is the optical lattice: a periodic 
scalar potential formed from standing waves of light [fl- 
optical lattices lead to simple bands with no topological 
character. In recent work[9], it has been shown that by 
using the optical fields to couple more than one internal 
atomic level, leading to effective gauge fields fl^, one can 
generate "optical fiux lattices" which break TRS and lead 
to bands with non-zero Chern number. Here we show 
that, by suitable optical dressing, one can form optical 
lattices that generate non-abelian gauge fields [sl. in 
a way that preserves TRS and leads to bands with non- 
trivial I2 invariant. Our approach is related to methods 



using optical fields to simulate the effects of spin-orbit 
coupling It is however very different to existing 

proposed cold-atom implementations of Z2 topological 
insulators [l2|. which imprint non-abelian gauge poten- 
tials on atoms in a deep optical lattice. Our key innova- 
tion is to approach the problem from the opposite, nearly 
free electron limit. This allows for large gaps, with size 
set by the recoil energy. Moreover, our proposal can be 
implemented with a relatively simple laser set-up. In ad- 
dition to its advantages to cold atom realizations, our 
nearly free electron approach brings a qualitatively new 
viewpoint to the field of topological insulators. 

We consider an atom with position r and momentum 
p and with N internal states which is described by the 
Hamiltonian 



Zm 



(1) 



where V has dimensions of energy. Ijy is the identity, 
while M{r) is a position-dependent N x N matrix acting 
on the internal states of the atom. To realize a Z2 topo- 
logical insulator requires N to be even, owing to a "spin" 
structure, and the Hamiltonian has to be invariant under 
time reversal, ^ = iCTj;i:[l|,[3]. (We denote by ai the Pauli 
matrices in spin space, and by K complex conjugation.) 
This requires that M = e-^Me. The smallest non-trivial 
case has = 4 with 



M 
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Here A,B,C are real numbers, and D 
real three-component vector. M can be also written as 



(3) 



where are the Pauli matrices in the block structure of 
Eq. (0). 

We propose to realize Hamiltonians of this form using 
four internal states of an atom, with the entries in Eq. ([2]) 
provided by optical fields Many implementations are 
possible. In principle the various transitions between the 
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internal states could be spectrally resolved, allowing the 
freedom separately to introduce coupling lasers for all 
transitions (or pairs of lasers for Raman coupling). Our 
results are very general, and could be applied to any such 
implementation. However, to make the discussion con- 
crete, we shall focus on a particular realization in which 
the coupling laser operates on a single frequency, and 
therefore has very limited freedom. 

The implementation we consider makes use of the 
properties of the groundstate (^5*0 — g) and long-lived 
excited state {^Pq = e) of ytterbium (Yb) or an alka- 
line earth atom. The usefulness of these two levels for 
optically induced gauge fields was pointed out by Ger- 
bier and Dalibard[13]: the long lifetime of the excited 
state reduces spontaneous emission; furthermore, there is 
a convenient "anti-magic" wavelength Aam at which the 
two states experience scalar optical potentials of oppo- 
site signs, allowing state-dependent potentials, ±Vam(^), 
to be implemented. We focus on Yb [3] which has 
nuclear spin / — 1/2. Then both g and e have two in- 
ternal states, and there are four states in total; these 
will correspond to the 4x4 structure in Eq. ([J). We 
consider the magnetic field to be sufficiently small that 
the Zeeman splitting is negligible, and all four e-g tran- 
sitions involve the same frequency wq = {E^ — Eg)/h. 
Combining the single-photon coupling, with electric field 



E = £e 



with the state dependent potential 



T4m, leads to the optical potential which in the rotating 
wave approxmation 15] is 
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where A = w — wo is the detuning. The form of the off- 
diagonal term follows from the Wigner-Eckart theorem, 
with dr the reduced dipole moment [16[. A comparison of 
Eq. dl]) with Eq. ([2]) shows that the optical coupling ^ 
describes a TRS situation provided all components of £ 
have the same phase, in which case £ can be chosen real. 

Two dimensional systems. We consider first cases in 
which the atoms are tightly confined in the z-direction, so 
their motion in the x—y plane is (quasi)-two-dimensional. 
A Z2 topological insulator can be formed by choosing the 
electric field, detuning and state-dependent potential in 
Eq. (jU such that 

dr£ — F (5, cos(r • Ki), cos(r • K2)) (5) 
^A + Km(r-) = VcOB[r-{Ki+K2)] (6) 

with Ki — (1,0, 0)k, H2 = (cos0,sin0, 0)k. The ampli- 
tudes are chosen to have a common energy scale F, which 
can be interpreted as a measure of the Rabi coupling. 
This energy scale should be compared to the characteris- 
tic kinetic energy scale, the recoil energy Er = ?i^K^/2m. 
The interspecies coupling ([SJ is formed from three (stand- 
ing) waves of linearly polarised light at the coupling 



frequency uj: two of equal amplitude with wavevectors 
in the 2D plane {ki for y-polarization and K2 for z- 
polarization) and one with wavevector normal to the 2D 
plane for a;-polarization with an amplitude smaller by a 
factor S. Since w ~ the magnitude of the in-plane 
wavevectors is k ~ 2tt/Xq with Ao = 578nm the wave- 
length of the e-g transition. The spatial dependence of 
Vam is set by a standing wave at the anti-magic wave- 
length Aam[13|, which creates a state-dependent poten- 
tial with -f K.2I — 47r/Aam- This fixes the angle 
e = 2arccos(±Ao/Aam). For Yb, Ao/Aam ^ 1/2 (to 
an accuracy of about 3%), so 6* ~ ±27r/3. For simplic- 
ity, in all following discussions we fix = 27r/3 and de- 
fine a = 47r/(-\/3«;). The optical coupling M then has 
the symmetry of a triangular lattice with lattice vectors 
ai = (V3/2, -l/2)a and aa = (0, l)a. 

The physics arising from this form of optical coupling 
is most clearly exposed by applying a unitary transfor- 
mation U = 2~^/^(l4 " «S3(T2)[17[ to the dimensionless 
coupling M 

M'^U^MU = ciEi-t-C2S2<5-3 + ci2E3-h(5S2CTi.(7) 

We have used the shorthand Ci = cosr • Ki and C12 = 
cos[r • {til + 1^2)]- For 5 — Q this matrix decouples 
into two 2x2 blocks, one for each eigenvalue of 0-3, as 
does the Hamiltonian since the kinetic energy ([T]) is di- 
agonal. Thus, the four-level system decouples into two 
two- level subsystems, each experiencing its own form of 
optical dressing, ciSi ± C2S2 + C12S3, differing only in 
the sign of the C2 term. This is precisely the optical 
dressing required to realize the "triangular" optical flux 
lattice described in Ref.[§] (See 18|). This optical flux 
lattice causes the atoms to experience a (periodic) effec- 
tive magnetic field with = 2 flux quanta in the unit 
cell. For V > 0.2i?fl the lowest energy band (which is 
twofold degenerate 0) is separated in energy from higher 
bands. This band is topologically non-trivial, having a 
Chern number of ±1, with sign set by the net sign of 
the coefficients ci,C2,ci2. Since this sign differs for the 
two subsystems, their lowest bands have equal and op- 
posite Chern numbers. (This is required by the fact that 
these subsystems are related by time-reversal.) Thus, 
for (5 = the optical dressing realizes the bandstruc- 
ture required to generate the quantum spin Hall (QSH) 
effect: there are two decoupled subsystems, the lowest 
energy bands of which have Chern numbers ±1. Each 
of these bands is filled when the 2D density of fermions 
(in that component) is equal to the magnetic flux den- 
sity, N^/{V^/2a?) = 4/(V3a2). When both bands are 
flUed, at a total fermion density of n2D = 8/-\/3a^, the 
bulk of the system is gapped. However, there will be a 
metallic surface state involving counterpropagating edge 
modes of opposite spin. This metallic surface state is the 
hallmark of the topological order that characterizes the 
QSH effect. 

The QSH effect is a special case of the Z2 topological 
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FIG. 1: Lowest energy bands for the 2D topological insu- 
lator with V = 0.5Er and 5 — 0.25 for wavevectors on the 
perimeter of the parallelogram (shown in the inset) with cor- 
ners r„m = (n/ti -|-mK2)/2 {n,m — 0,1). Each band is twofold 
degenerate, with fourfold degeneracies at r„m 7^ 0. The four 
bands with E/Er < —0.1 have non-trivial Z2 invariant, as 
follows from the eigenstate parities (±) at r„m. 

insulator. The Z2 topological insulator, in general, does 
not allow separation into two decoupled systems, but re- 
tains counterpropagating edge states that are protected 
only by TRS. A nonzero value of 6 introduces a coupling 
between the two spin-states in a manner that preserves 
TRS, in direct analogy to the effects of spin-orbit cou- 
pling on the electronic structure of solids. In FiglT] we 
show the few lowest energy bands for a small nonzero 
value of 6. (The bands were calculated by numerical diag- 
onalization in the plane wave basis.) For S — all bands 
are fourfold degenerate. For 5 7^ this degeneracy splits 
into two doubly degenerate branches at most wavevec- 
tors. However, the fourfold degeneracy remains at three 
nonequivalent symmetry points in fc-space: Ki/2, 
and (ki -I- H2)/2. (The reciprocal lattice vectors are sim- 
ply given hy Ki.) 

Due to the coupling between the two subsystems, the 
Chern number of a subsystem is no longer well defined, 
but the bands still carry a topological character since the 
6 = case is adiabatically connected to the QSH limit. 
To demonstrate this, one has to use a more general Z2 
invariant [loj. This invariant takes a particularly simple 
form for our system owing to its symmetry under spatial 
inversion, P : r — )■ —r. Inversion symmetric insulators, 
as shown by Fu and Kanefiol [20j. are Z2 nontrivial if 

n n (8) 

n,m— 0,1 QGfillcd 

where are the inversion eigenvalues corresponding 
to the a-th Kramers doublet of bands at the momenta 
r„m = {nKi + mK2)/2 {n, m = 0, 1). Evaluating Eq. ^ 
for our system we find that it is indeed in a topological 
phase: the points with fourfold degeneracies contribute 
with (—1) while fc = contributes with unity, as in- 



FIG. 2: Lowest energy bands for the model Eq. ((9]), for a set 
of wavevectors k = kiKi with k2 and kz uniformly spaced 
across the Brillouin zone. The parameters are V = Q.QEr, 
5 = 1.0, ^i = -0.2. The four bands with E/Er < -0.9 have 
non-trivial Z2 invariant which, together with the clear band 
gap, demonstrate the topological insulator phase. The right 
hand panel shows the density of states. 

dicated in Fig [TJ Although we have shown bands for 
the special case of a triangular lattice, our approach is 
very general: starting from an optical flux lattice of any 
symmetry our method will lead to a 2D Z2 topological 
insulator. 

Three dimensional systems. Our approach has the re- 
markable feature that it generalizes naturally to 3D topo- 
logical insulators. One simply needs to substitute 

C12 ^ C12 + S{fi + C13 + C23) (5 ^ (5cos(/t3 • r), (9) 

where our notation follows that of Eq. ([7]), and we use 
/«3 = (0,0, That this system is a good candidate for 
a topological insulator can be checked using the parities 
as in Eq. ([5]) , but the product is now over eight momenta 
r„m/ = {nKi+mK2+lK3) /2 (n, m, Z = 0, l)tl9|. Evaluat- 
ing the product shows that the parities are those of a Z2 
topological phase for a large range of parameters. Impor- 
tantly, this range includes systems which are insulators 
as well. This is demonstrated in Fig. [5J where we show 
the few lowest energy bands of the model for a point in 
the regime where the parity product is negative. A clear 
bandgap separates the lowest four bands from the higher 
bands, showing that the system is an insulator. This ex- 
tension to 3D entails additional geometrical constraints 
on the wavevectors of the coupling and state-dependent 
lasers. For example, for Ki _L Kj the state-dependent po- 
tential Cij has to have periodicity Ao/a/2. For Yb, this 
differs from that achieved by the anti-magic wavelength 
Aam/2 ~ Ao, so in this case there will be an additional 
state independent potential We have checked that 
this does not affect the appearance of the Z2 topological 
insulator phase for the lowest 4 bands. 

While the robust and simple realization of a 3D topo- 
logical insulator is a key result, the choice ^ might seem 
fortuitous. The remaining part of the Letter introduces 
a constructive approach for finding 3D systems with Z2 
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nontrivial parities. A hint for this is given by the 2D 
case, if one notes that the Tnm with negative parities are 
precisely the points in momentum space with fourfold de- 
generacies. This can be understood using the following 
symmetry considerations. The Hamiltonian is invariant 
under translations by the lattice vectors a^. The asso- 
ciated reciprocal lattice vectors are simply Kj, and de- 
fine the Brillouin zone for the conserved momentum fc. 
Using the Bloch decomposition i/'fc — exp(zk • r)ufc(r), 
the eigenvalue problem at each fc point follows from the 
Hamiltonian 

H^^^^^lN + VM'{r) (10) 
acting on lattice periodic functions. has half- 



translation symmetries Ti — Ti2T^ 



ai/2 



and To — 



Yiia3Ta^/2. (Here Tv is the operator of translation 
with the vector v.) At the special momenta r„m, 
Hr„m is invariant under both generalized time-reversal 
and inversion operators: 9nm = exp(— i2r„m • r)9, and 
Pnrn = exp(— i2r„m • r)P. Thc two symmetries com- 
mute, hence ur„„ and its Kramers partner 0nmUr„,„ 
share the same inversion eigenvalue ^nm- From the fact 
that = 1 on lattice periodic functions, it follows that 
exp(-ir„m • ai)PnmTi. This means that at Tio 
the application of Ti, at Fqi of T2 and at Fn of either 
Ti reverses the inversion eigenvalue: this requires, firstly, 
that the energies are fourfold degenerate, and, secondly, 
that for each energy the inversion eigenvalues come in op- 
posite pairs 21| . This explains why the degenerate points 
contributed with (—1) in Eq. ([8]). 

These considerations allow one to construct 3D con- 
figurations with Z2 nontrivial parities: consider a model 
which has a reflection symmetry ai -H- 02 and a half- 
translation symmetry by d = i (01 + 02 + 03). (03 is the 
third primitive lattice vector.) The reflection symmetry 
results in ^nmi = Cmni, while the half translation sym- 
metry means that at Fooi and Fm the inversion eigen- 
values come in opposite pairs for each energy, making 
the energies fourfold degenerate. Thus, for any 2N dou- 

1, so the net parity is 



blets we have n^^i Unm Cimi 



a = l 1 Inm >nmO' 



Assume that the model can be written 



(pert) 



as Hf, = 42^) + (^±^ r 

k 2m fc 

a two dimensional topological insulator with 2N filled 



If H^^^^ realizes 



Kramers pair of bands, then Ha^i Dnm^ 
a weak Hjf"^\ This means that the net parity for the 
lowest 2A'' Kramers pairs must be negative: thc model 
is guaranteed to have the parities of a 3D Z2 topological 
insulator. The choice Eq. ([SJ is one of the simplest ones 
consistent with this strategy. 

In summary, we have described how simple forms of 
optical coupling can lead to low-energy bands that have 
non-trivial Z2 topological invariant, in both 2D and 3D. 
Weakly interacting fermions filling these bands will form 
incompressible band insulators, and exhibit the expected 



f(a) _ 
nmO 



-1 for 



features of these topological phases: notably the metal- 
lic surface states. The properties of these surface states 
could be probed by studying thc (spin-resolved) collec- 
tive modes of the atomic gas. Since we work close to the 
nearly free electron limit, the required temperature scale 
is set by the recoil energy, which is of the order of 0.1 mK 
in thc case of -'^''-'^Yb, easily accessible in experiments. 

The realization of these topological bandstructures in 
ultracold atomic gases will allow the study of interest- 
ing properties that cannot readily be explored in solid 
state systems. An important issue concerns the effects 
of (attractive or repulsive) interactions (22] . Weak inter- 
actions can affect the properties of the surface modes, 
while strong interactions may drive the system into 
strongly-correlated topological phases, either with TRS 
preserved [23}. or with broken TRS in which case they 
are related to fractional quantum Hall states. Further- 
more, cold atom implementations will allow studies of the 
effects of tuning the nature of the bands, inducing tran- 
sitions in their topology. Introducing a phase difference 
between the electric fields in (j4]) breaks TRS; this could 
be used to study quantum phase transitions out of the 
Z2 topological insulator phase, into (for example) integer 
quantum Hall states. 
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SUPPLEMENTARY MATERIAL 

In this Supplementary Material, we briefly summarize 
some basic results from Ref. [l^ about evaluating the Z2 
invariant in inversion symmetric systems. We repeat 
them here to help integrate the key formula, Eq. 
below, into the literature on cold atomic gases. A more 
elaborate discussion including detailed derivations can 
be found in Refs. [H, 111. 

A time-reversal invariant band-insulator is character- 
ized by a Z2 invariant v, with 1/ = and i' — 1 denoting 
trivial and non-trivial insulator respectively. The value 
of this invariant follows from the properties of the Bloch 
wavefunctions of the occupied bands. For the generic case 
(with no symmetries other than time- reversal), there are 
several equivalent ways for its evaluation. However, all 
methods require a careful gauge choice for the wavefunc- 
tions across the Brillouin zone. Finding such a gauge is 
often difficult. 

Shortcuts are possible in the presence of additional 
symmetries. For example, for 2D systems where the z- 
component of the spin is conserved, evaluating v reduces 
to the calculation of the Chern number C from the wave- 
functions of one of the spin components. The result is 
simply 

i^ = Cmod2. (11) 

This result is used to establish the quantum spin Hall 
phase of Eq. ([7]) of the main text for i5 = 0. 

In the main text, we extensively use another shortcut, 
which is available if the system has inversion symmetry. 
This leads to a radical simplification in both 2D and 3D. 
Before citing the result of Ref. [Til we briefly summarize 
some basic facts about the eigenstates of inversion sym- 
metric, time-reversal invariant, lattice periodic Hamilto- 
nians. Inversion symmetry means that the Hamiltonian 
H{r) is invariant under P : r ^ — r. The Bloch Hamil- 



tonian is 

i?k(r-) - e-*'='^i^(r)e''^^ (12) 
and it acts on lattice periodic functions (Bloch functions) . 
We denote the reciprocal lattice vectors by Kj. At a 
generic fc, unlike H{r), the Bloch Hamiltonian is not 
invariant under P nor under time-reversal 9, since both 
map k — >■ —k. The momenta r{„^.} = njHj /2 (rij — 
0,1), however, are special because T^nj} and — r{„^.j. 
are equivalent wavevectors. This means that 77(r)r{„ .} 
is invariant under both generalized time-reversal and 
inversion operators: ^{nj} — exp(—z2r{„^ }r)6', and 
{" } ~ 6xp(— i2r{-„ .}r)P, where the exponentials take 
Hir)-ri^.y back to k(r)r{„^.} [cf. Eq. Note that 

the operators d{n } ^^id P{n } rnap lattice periodic func- 
tions to lattice periodic functions. 

Because of the generalized inversion symmetry, the 
eigenstates u{r)r^^.y of H{r)r^„.^ can be labeled by in- 
version eigenvalues £,{nj} — ilj 

A",}"(^)r{„^} =^{„^}u(r)r{„^}. (13) 

It is straightforward to show that the corresponding 
eigenstates il'{r)r^„.y = exp(jr{„^.}r)M(r)r{„ .j of the 

Hamiltonian H{r) are eigenstates of P (and not P{nj}) 
with the same eigenvalue ^{n^}- 

The symmetry under ^{,1^} means that the energies 
at r{„^.} (i.e., the eigenvalues of /f(r)r{„.}) come in 
Kramers degenerate pairs. The Kramers pair of a 
Bloch eigenstate zi(r)r{„ .} is 6'{nj}u(T')r,„ j • Because 

the symmetries 6{nj} a-nd P{nj} commute, u{r)^^^.y and 
0|„j}u(r)r{„ .} share the same inversion eigenvalue. The 
a-th Kramers doublet at r{nj } can thus be characterized 
by a single inversion eigenvalue 

We are now ready to state the result of Ref. 119|: Fu 
and Kane showed that the inversion eigenvalues corre- 
sponding to the filled bands straightforwardly determine 
whether an insulator is in a topological phase. The sys- 
tem is a topological insulator (v = 1) if 

n n ^12}=-!' (14) 

{rij} aefillcd 

where a runs over the Kramers doublets of the filled 
bands. Note that the invariant is even simpler than 
the Chern number, since it does not require the integra- 
tion of a (possibly complicated) function over the Bril- 
louin zone; its sole inputs are the parity eigenvalues at a 
few discrete points in momentum space. 



